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Atomic masses of the neutron-rich isotopes "^"""^''Cd, "^^'"'In, """"^'Sn, """""Sb, and 
i32-i40rpg iiave been measured with high precision (10 ppb) using the Penning trap mass spec- 

^^^Sn, i3«Sb, and 139,140^^ 



trometer JYFLTRAP. Among these, the masses of four r-process nuclei Sn, Sb, 
were measured for the first time. An empirical neutron pairing gap expressed as the odd-even stag- 
gering of isotopic masses shows a strong quenching across N—82 for Sn, with the ^-dependence that 
is unexplainable by the current theoretical models. 



PACS numbers: 21.10.Dr, 21.60.-n, 27.60.+J 

The doubly magic ^'^^Sn nucleus has been probed 
intensively by nuclear spectroscopy over the last two 
decades. It has been found to exhibit features of ex- 
ceptional purity for its single particle structure [l|, y]- 
This provides an ideal starting point for exploring de- 
tailed evolution of nuclear structure of more neutron-rich 
nuclei beyond the N=82 closed shell in the vicinity of Sn. 
Only a few experimental and theoretical attempts along 
these lines have been performed recently. No experimen- 
tal data exist for excited states or masses for nuclides 
below Sn with A^ > 82. The experimental situation is 
slightly better for Z > 50 isotopes of Sb and Te because 
of their easier access. 

Recent data on the B(E2) transition strengths for 
i32Te, 134^6 and ^^^Te isotopes [3] and their interpre- 
tation using a quadrupole-plus-pairing Hamiltonian and 
Quasiparticle Random Phase Approximation (QRPA) [4] 
suggested the need for reduced neutron pairing to explain 
the observed anomalous asymmetry in the B(E2) values 
across the N=82 neutron shell. This behaviour _was not 
observed in standard shell model calculations [SJ. An- 
other shell model calculation of the binding energies of 
heavy Sn isotopes with A > 133 [5| suggested the im- 
portance of pairing correlations and the strength of the 
pairing interaction in general for weakly bound nuclei. 
Therefore, it would be necessary tq^ probe the evolution 
of odd-even staggering of masses [6| around the A^ = 82 
neutron shell to learn about the magnitude of pairing and 
its variation as a function of Z and N beyond ^^^Sn. 

High precision of present-day ion-trap mass spectrom- 
etry combined with high sensitivity [7| can provide the 
needed information on mass differences such as one- and 
two-nucleon separation energies, shell gaps and empirical 
pairing energies. For example, the masses of neutron- 
rich Sn and Xe isotopes were recently measured up to 
^^''Sn and ^''^Xe with a Penning trap mass spectrome- 
ter ISOLTRAP at the CERN ISOLDE facility f8','9l|. In 
this Letter we wish to present new data of high-precision 
mass measurements of neutron-rich Cd, In, Sn, Sb, and 
Te isotopes across the A^=82 neutron shell by using the 



JYFLTRAP Penning trap. These nuclides are also of 
interest for nuclear astrophysics models of element syn- 
thesis, in particular, to explain the large r-process abun- 
dance peak at A=130 [10|, see Fig. [TJ In more general 
context, a vast body of nuclear data on neutron-rich iso- 
topes is needed for r-process nucleosynthesis predictions. 
Such data include masses, single particle spectra, pairing 
characteristics as well as decay properties and reaction 
rates. In all of these the binding energies or masses of 
ground, isomeric and excited states play key roles |lO| . 
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FIG. 1. Neutron-rich isotopes with Tz > 13 whose masses have 
been determined by Penning trap (PT) or Q^-measurements. 
Letter r denotes r-process nuclei according to Ref. \W^ . 



The measurements were performed using the 
JYFLTRAP Penning trap mass spectrometer |ll| 
which is connected to the Ion Guide Isotope Separator 
On-Line (IGISOL) mass separator [12|. The ions of in- 
terest were produced in proton-induced fission reactions 
by bombarding a natural uranium target with a proton 
beam of 25 MeV energy. A thorium target was used in 
the case of ^^^In and isotopes of Sb. 

Fission products stopped in a helium-filled gas cell at 
a pressure of about 200 mbar as singly-charged ions were 
transported out of the gas cell, accelerated to 30 keV 
energy, and mass separated. A gas-filled radio frequency 
quadrupole cooler and buncher prepared the ions for the 



Penning trap setup. 

In a Penning trap an ion has three different eigenmo- 
tions: axial motion (v^) and two radial motions, mag- 
netron (i^_) and modified cyclotron (i^+) motion. The 
frequencies of the radial motions sum in first order to 
the cyclotron frequency i/c = ^^B. Here q and m are 
the charge and mass of the ion, and B is the magnetic 
field. 

JYFLTRAP consists of two cylindrical Penning traps, 
the purification and precision trap, which are located in- 
side a 7-T superconducting magnet. Both traps were 
used to purify the samples, first with the sideband cool- 
ing technique [l3| for isobaric separation and then, if 
prompted by the presence of isomers or other contam- 
inants, with the Ramsey cleaning technique |l4l |. After 



a cleaned sample was trapped in the precision trap, the 
time-of-flight ion-cyclotron resonance (TOF-ICR) tech- 
nique |15| was applied in order to determine the reso- 
nance frequency. The Ramsey method of separated os- 
cillatory fields jig , which makes the sidebands more pro- 
nounced and narrower, was used in the majority of mea- 
surements. The experimental setup and the measure- 
ment technique are described in more detail in Ref. |17J . 
A sample resonance is shown in Fig. [2j When the cy- 
clotron frequencies of an unknown species (wmeas) and a 
well-known reference ion (mref) are known, and both are 
singly charged ions, the atomic mass can be determined: 

Wmeas = — ^^^^— (™rof -Tie) +me, (1) 

^c,mcas 

where rUe is the mass of the electron. The data analysis 
procedure was almost identical compared to [18| . In this 
work, the systematic uncertainties due to the magnetic 
field fluctuations were minimised by applying so-called 
interleaved frequency scanning described in Ref. [19|. 




-6-4-202468 
Excitation frequency - 796410.46 [Hz] 



FIG. 2. (Color online) A time-of-flight (TOP) resonance from 
the JYFLTRAP setup for "^Sn+. A two- fringe Ramsey pat- 
tern of 25-350-25 ms (on-off-on) was used in this case. 

The results are given in Table HI Excitation times be- 
tween 100 ms and 800 ms were used depending on the 
half-life of the isotope. In the case of Ramsey excitations, 
excitation patterns with two 25 ms pulses separated by a 
waiting time, were used. Table U lists only the masses of 



the ground states which are relevant for the further dis- 
cussion of the results. A paper containing the isomeric 
data will be submitted separately. The new data agrees 
with the earlier ion-trap measurements and presents sig- 
nificant improvement in accuracy for all Sn, Sb and Te 
isotopes beyond A^=82. 



TABLE I. Cyclotron frequency ratios f and ground state mass 
excess values in keV based on this work. Masses of reference 
Xe isotopes in column 1 for given A are from [20l - [22 |. Results 
from other direct mass measurements are from ISOLTRAP @, 
|23i[2JI or an experimental storage ring (ESR) at the fragment 
separator (FRS) [25^ 
Evaluation 2003 (AME2003) H is given 



Otherwise, the value from Atomic Mass 
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and filling approximations [271 128| were used to treat odd 
nuclei. 
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FIG. 3. (Color online) Experimental odd-even mass stag- 
gering for Sn, Te, and Xe isotope chains. For clarity, odd- 
A^ points have been connected with lines and the error bars 
have been omitted when 5(A'^') < 10 keV. The open symbols 
present the experimental values around the shell-gap prior to 
this work. The data points for N=82 are beyond the scale. 

The high precision of Penning-trap mass measurements 
enables for the first time a critical evaluation of odd-even 
staggering of binding energies and related empirical pair- 
ing gaps across the N = 82 shell gap. The most simple 
example is the three-point odd-even-mass-staggering for- 
mula [a| 

A^'^\n) = {-l)^[E{N + l)-2E{N) + E{N-l)]/2, (2) 

where N is the number of nucleons (neutrons or protons) 
and E the binding energy. The A^'^^ staggering mostly 
depends on the intensity of pairing correlations in nu- 
clei, but, as we discuss below, it is also affected by the 
polarisation effects. 

Figure [3] shows the experimental neutron A*^^^ stagger- 
ing for Sn, Te, and Xe isotopes crossing the N=82 shell 
closure. The difference between the values at N—81 and 
83 shows a large asymmetry for Sn but a much smaller 
one for Te and Xe. This indicates a considerably stronger 
quenching in pairing gap for Sn than for Te and Xe sug- 
gesting the importance of core polarisation effects. A 
similar asymmetry observed for B(E2) values of n-rich 
Te isotopes was also traced to reduced neutron pairing 
above the N=82 shell closure 0]. 

In order to probe this question theoretically, we per- 
formed self-consistent calculations of varying complexity, 
by using the Sly4 ^261] energy density functional and con- 
tact pairing force. To make a meaningful comparison, 
in each variant of the calculation, the pairing strength 
(equal for neutrons and protons) was adjusted so as to 
give the average pairing gap in ^^"Sn equal to 1.245 MeV. 
The pairing channel was described within the Hartree- 
Fock-Bogoliubov (HFB) approximation and the blocking 
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FIG. 4. (Color online) Same as in Fig.Obut for the calculated 
(spherical HFB) odd-even mass staggering in (a) Sn, Te, and 
Xe isotope chains and mixed pairing force, and (b) Sn isotope 
chain and surface, mixed, and volume pairing forces. 

First, to provide a baseline for further analyses, in 
Fig. [4] we show the neutron A'^'^^ staggering calculated 
within the spherical approximation [29]. Such an ap- 
proximation allows us to look at pure effects of pairing 
correlations, whereby the binding energies of odd iso- 
topes are decreased due to one pair being broken. Oth- 
erwise, in this approximation, effects due to deformation 
polarisations exerted by valence particles are completely 
switched off. From Fig. Hb we see that for the volume 
(vol) or mixed (mix) pairing forces [30J, the experimen- 
tal decrease of A^^^ when crossing the N=82 gap in Sn 
is very well reproduced, whereas the data exclude the 
pure surface pairing force. Such pairing decrease is due 
to a lower level density above the N=82 gap [4| (the sur- 
face pairing force is unable to discriminate between the 
two surface-type orbitals vhn^ and vir/2 located below 
and above the N=82 shell gap, respectively). Spherical 
calculations miss the experimental values in Te and Xe 
isotopes (Fig. [4^) , which indicates that the polarisation 
effects must here be taken into account explicitly. 

To better illustrate the trends in Z > 50 nuclei, in 
Fig. [5] we show the experimental neutron A'^'^^ stagger- 
ing in Sn, Te, and Xe isotopes along the iV=81 and 83 
lines, compared with theoretical results. In experiment 
(Fig. [5K), the difference between the N=81 and 83 iso- 
tones smoothly decreases from about 0.5 MeV in Sn to 
almost zero in Xe, whereas spherical results (discussed 
above and repeated in Fig. [SJa) show no such a decrease 
at all. 

To analyse the effects of polarisations induced by de- 
formation, we have performed two additional HFB cal- 
culations. First, by using the code HFODD (v2.51i) 
|3l| . we allowed valence particles or holes to induce self- 
consistently deformed shapes. Then, the even isotones, 
N — 80, 82, and 84, turn out to be spherical anyhow, 
namely, neither the closed-core N — 82 systems become 
deformed, nor the paired two neutron particles or holes 
induce non-zero deformations. Only the polarisation ef- 
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FIG. 5. (Color online) Odd-even staggering for the A'^ = 81 
and 83 isotones as measured in experiment (a) and estimated 
within the spherical (b), deformed (c), and deformed particle- 
number-conserving (d) self-consistent calculations. 



fects exerted by unpaired odd neutron particles or holes 
are strong enough to induces non-zero deformations in 
the odd isotones, A^=81 and 83. 

In this way, the odd-particle polarisations lead to lower 
values of the A^'^^ staggering solely through increased 
binding energies of odd isotones. Since such polarisa- 
tion increases with adding protons, the obtained values 
of A^^^ decrease with Z, as illustrated in Fig. [St. How- 
ever, this trend does not depend on whether deformation 
is induced by odd particles or holes; therefore, on both 
sides of the N=82 shell gap we obtain identical decrease 
with Z , at variance with experiment. 

To test if the above results may be affected by the 
particle-number non-conservation, inherent to the HFB 
theory, we have repeated all calculations by using the 
HFBTHO code [13, IHl. Here, within the Lipkin- 
Nogami method, we were able to include the approximate 
particle-number projection after variation. Moreover, 
the obtained solutions were next exactly projected on 
good particle numbers. Such a projected Lipkin-Nogami 
(PLN) method was extensively tested f33] and proved 
to be very efhcient in describing pairing correlations in 
near-closed-shell systems. 

The obtained results are shown in Fig.[5ji. We see that 
the general pattern is not qualitatively changed. The dy- 
namic pairing correlations, which are now nonzero even 
in the N=82 isotones, lead to larger values of the A'^' 
staggering, which in the iV=81 isotones perfectly well re- 
produce the experimental trend. However, in the N=83 
isotones, the disagreement with data remains a puzzle. 
Indeed, the asymmetry of the trend of the staggering, 
measured below and above the A^=82 gap, points to spe- 
cific effects related to orbitals occupied beyond N—82 or 
to the their weak binding, which are not captured by the 
current state-of-the-art theoretical approaches. 

This work has been supported by the Academy of Fin- 
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